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The fundamental quantum dynamics of two interacting oscillator systems are studied in two 
different scenarios. In one case, both oscillators are assumed to be linear, whereas in the second 
case, one oscillator is linear and the other is a non-linear, angular-momentum oscillator; the second 
case is, of course, more complex in terms of energy transfer and dynamics. These two scenarios 
have been the subject of much interest over the years, especially in developing an understanding 
of modern concepts in quantum optics and quantum electronics. In this work, however, these two 
scenarios are utilized to consider and discuss the salient features of quantum behaviors resulting from 
the interactive nature of the two oscillators, i.e., coherence, entanglement, spontaneous emission, 
etc., and to apply a measure of entanglement in analyzing the nature of the interacting systems. 

The Heisenberg equation for both coupled oscillator scenarios are developed in terms of the rele- 
vant reduced kinematics operator variables and parameterized commutator relations. For the second 
scenario, by setting the relevant commutator relations to one or zero, respectively, the Heisenberg 
equations are able to describe the full quantum or classical motion of the interaction system, thus 
allowing us to discern the differences between the fully quantum and fully classical dynamical pic- 
ture. 

For the coupled linear and angular-momentum oscillator system in the fully quantum-mechanical 
description, we consider special examples of two, three, four-level angular momentum systems, 
demonstrating the explicit appearances of entanglement. We also show that this entanglement 
persists even as the coupled angular momentum oscillator is taken to the limit of a large number 
of levels, a limit which would go over to the classical picture for an uncoupled angular momentum 
oscillator. 
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1. INTRODUCTION 

In this paper, we study the quantum dynamics of two interacting oscillator systems in two different, but very well- 
known scenarios. In case one, both oscillators are assumed to be linear, and in case two, one oscillator is assumed to 
be linear, while the other oscillator is assumed to be non-linear (angular momentum oscillator). We study Heisenberg 
equations and the time evolution of the wave-functions \ip(t)) = e - lHt / h \yb{0)) for both cases to discuss the quantum 
interactive behaviors between the two oscillators therein, e.g., entanglement, spontaneous emission (in case two), and 
quantum versus classical differences (in case two). Entanglement, an implicit non-classical ingredient of interacting 
quantum systems, is the key to enabling quantum computation, which allows for the solution of certain classes of 
problems more efficiently than the classical counterpart (see, e.g., Refs. HQ). For case two, by setting the relevant 
commutator relation in each oscillator to one or zero, we utilize an interesting quantum-classical schemed to clearly 
characterize the quantum features in the system. 

In order to investigate the time evolution of the wave-function and the notion of entanglement in case one, we 
adopt the coupled-boson representation^ which was discussed by Schwinger in terms of the correspondence between 
two one-dimensional linear oscillators and an angular momentum oscillator. Based on this representation, we will 
easily obtain exact energy eigenvalues and eigenstates of case one. The coupled-boson representation has been used in 
context with similar types of problems to construct an ideal model for quantum-mechanical interference^ and to study 
the current operator in two coupled linear oscillator (the rate of exchange of occupation number between oscillators), 6 

We adopt here 1-Tr {p {v) ) 2 as an entanglement measure for the total system in a pure stated where p^ v ' is a reduced 
density matrix of individual oscillator v. We will show that this measure increases with the occupation number n for 
a given initial state for case one, and with the maximum quantum number j of the angular momentum oscillator for 
case two, respectively; this increase with occupation number results from considering a dissipationless and coherent 
process, whereas we know that entanglement will diminish when loss is taken into account, e.g., resulting from the 
(thermal) interaction between the system and its environment.— For case two, we will discuss entanglement for the 
cases where the angular momentum oscillator is in the states j = \, 1, § explicitly. 

In section [21 we study the interaction between two linear oscillators. The Heisenberg equations of motion are 
analyzed in the coupled boson representation, leading to the straightforward diagonalization of the coupled oscillator 
Hamiltonian; from this diagonalized Hamiltonian, the exact eigenstates and eigenvalues of the coupled oscillator 



2 



Hamiltonian are obtained, thus allowing for the calculation of the reduced density matrix, and the subsequent study 
of the entanglement under the two initial conditions of simple product state (a disentangled state) and exact eigenstate 
for the coupled oscillator system (an entangled state). 

In section we study the interaction between a linear and a non-linear, angular momentum oscillator. We make 
use of the seminal work of Ref. to analyze the quantum-mechanical and the classical dynamics of this interaction, 
and to connect the quantum manifestation of spontaneous emission in this system with the notion of entanglement. 
Further, the total system Hamiltonian is diagonalized for j — ^, 1, | using rotational generating operators of the 
group SU(n)£ From this diagonalization, the reduced density matrix, and thus the measure of entanglement, is 
calculated for simple product (disentangled) initial states. Unlike the linear-linear oscillator system, this coupled 
system displays periodicity in the measure of entanglement for j = \ and j = 1, whereas for j = | and beyond, 
the measure of entanglement is aperiodic; this aperiodicity is shown to be a result of the complex structure of the 
diagonalized multi-level Hamiltonian. 



2. INTERACTION OF TWO LINEAR OSCILLATORS 
2.1. Quantum behavior in the Heisenberg picture 



Let us begin with case one. The coupled oscillator system under investigation is described in the rotating wave 
approximation by the Hamiltonian 

H = huj\ ( a\ di + ) 4- hui2 (4 0.2 + ) + h n ( a\ 0,2 + Qi 4 ) i (1) 



where n denotes a coupling strength, and we assume that uj-y > uj 2 . Schwinger found that the quantum-mechanical 
angular momentum can be obtained by using creation and annihilation operators of two one-dimensional harmonic 
oscillators in terms of Jl oc a' ® a. This is explicitly given as follows: 4 

J x = \ (a\ ® a 2 + &i <g> 4) , Jy = i (a\ ® a 2 - o x <8> 4) > Jz = | ("1 " "2) , J = f (™i + n 2 ) , (2) 

where h v — b%a v with v = 1,2 is obviously an occupation number operator of each linear oscillator v. From 
[fip , a\\ = 1 5^1, the operators Ji, where I = x,y, z, satisfy the usual angular-momentum commutation relations. The 
operator J yields the quantum number of angular momentum, j = 0, |, 1, 1, • • • . In terms of the angular momentum 
operators in J2]) the Hamiltonian in {1} is rewritten as 

H = Auj J z + 2kJ x + (uji +oj 2 ) J + |fi {0J1 +uj 2 ) , (3) 

where Auj = u)\ — uj 2 > 0. Since [jl, J] = 0, we easily obtain a constant of motion J with j = \{ni + n 2 ) : 



ihj — 



0, (4) 



which means that the total oscillator occupation number is conserved. Also, we find that J x — ~AuiJ y and so J x is 
a constant of motion for Aui — 0. 
With the aid of the identity, 

e w v e -*0 = v + V } + f [L>, [L>, V]] + f [U, [U, [U, V}}} + ■■■, (5) 

we evaluate e %lJy ^/ h H e~ llJy ^/ h , where < 7 < -|, to arrive at the diagonalized Hamiltonian of eq. |(3]) as 

H d = e ^^' h H e -niy(0)/ft = | % ( Wl + w 2 ) + {Auj ■ (cos 7 ) + 2k ■ (sin 7)} J z (0) + + uj 2 ) J , (6) 

where tan 7 = ^ (note that in case of Auj = 0, we have 7 = § for any non-zero k). From eqs. ([2]), |(6]) we now obtain 
two renormalized uncoupled harmonic oscillators with the renormalized energy splitting for each oscillator, 

H d = huj[ (A\ A x + i) + huj' 2 (a\A 2 + I) (7) 
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where = \{ux+<jj 2 ) + \u) ,u' 2 = \{ui+w 2 )-\u>. Here, d 2 = (2k) 2 + {Alu) 2 , and A v (t) = e-^WA a v {t) e i ^^/ n 
explicitly reads 




where the rotation matrix is 

A/ ,n ( cost) -siwd \ , , 

R W = [sin$ cos ) ■ ® 

Thus, from the Schwinger approach, we realize that the coupling term in eq. (U) can be equivalently represented as 
the J x component of angular momentum, and therefore, H in eq. (1} can be diagonalized by a unitary rotational 
operator with respect to the angular momentum operator J y {Q) through 7.— Then, using eqs. ([7]), iJH) we easily obtain 
the original annihilation operators as 

fi„(t) = {(cos^) 2 e-*"'"* + (sin^) 2 e~H f } a„(0) + \ (sin 7) (e^ 4 - e- <w =*) a M (0) , (10) 

where v,H= 1,2 but v ^ \i. Here, neglecting the "non-local" observable a M (0), we obviously have [&!/(£), aj,(i)] 9^ 1 at 
i 7^ 0. Now, we easily acquire the exact eigenstates and eigenvalues of the Hamiltonian in (U) as 

\n u n 2 ) H = e-^W/* |m) |n 2 ) ; £„ 1)Tl2 = § (m + n 2 + 1) (wi + c^ 2 ) + § (m - n 2 ) w . (11) 

We see here that for 7 = 5 (i.e., Aw = 0), the eigenstates |ni,n 2 ) ff will also be eigenstates of j x (0) = J x {t)- 



2.2. Appearance of the entanglement in the Schrodinger picture 

Let us now study the entanglement in the time evolution of the two coupled linear oscillators, \ip(t)) = e~ lHt ^ h \ip(0)). 
To this end, we adopt,— as a measure of entanglement, 

M m)) = 1 -7>[pW(t)] , (12) 

where V \p- v '\ = Tr„ (p^) 2 is the purity measure of the state of each oscillator v. Here, it clearly holds that 
■pf/S*- 1 - 1 ] — V[p^]i^ As can also easily be verified, V[p^ v '\ takes on a maximum value of 1 if the oscillator v is in a pure 
state; it then follows that M\^) — for a simple product (disentangled) state, \tp) = \tfj^)\(/)^) . On the other hand, 
V[p^] takes on its minimum value -p, where M is a number of accessible orthogonal states {\pv);p = 1, 2, • • • , M} 
of each oscillator v, if the oscillator v is in the maximally mixed state given by p^ = jj 5 k i] then, the corresponding 
total wave-function is a maximally entangled state in form of -7= (|li) 1 12) + |2x) |2g) + ••• + |Mi)|M2)), and 

has its maximum value . We will explicitly determine M\^( t )) in its time evolution with mathematically 
manageable initial states |^>(0)) given below. 

We obtain the reduced density matrix, p( v \t), from the system state evolution, \ip(t)}, by constructing the system 
density matrix pit) = \ip(t)) (ip(t)\ from the known state of the system, and then taking the trace of p(t) with respect 
to all subsystem matrix elements within the system but those of subsystem v. For our oscillator system under study, 
the state evolution is given by 

00 

|V»(*)> = c «i^ K>"2) ff e"^"!-"*', (13) 

m ,n2—0 

where c ni] n 2 — (^(0)| nx^n%) H . Therefore, in constructing p = \tp) from (fl3|) . and taking the trace of p with 
respect to oscillator 2, we formally acquire the reduced density matrix for oscillator 1 as 

oo 

pi)\t) = Tr 2 p(t) = ^{p(C ;ls . (14) 
where {p(i)} fes . ;s = (»2| (&i|p(£) \h) \s 2 ). Using \ip(t)) in eq. HH), we explicitly obtain 
Pu(t) = EE c — ^ • c k,n' 2 ■ e-^ E ^- E ^ ]t . (s 2 \ (fell e~^ n K) |n 2 ) • (n' 2 \ (n[\ e^' h \h) \s 2 ) , (15) 

"1,12 n' lt n' 2 s 
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where J y = J y {0), and \ni,n2)H from eq. ifTTj) has been used in eq. i|13p to obtain the explicit result. 

Now, in transforming eq. lfl5|) to the coupled-boson representation discussed in eq. (J2j) , we relabel |m)|ri2) and 
l n 'i)l n 2) m terms of the eigenstates of j z , namely, \j; to) and \f; m'), respectively, where j — |(ni+rt 2 ), m — |(ni— n 2 ) 
and j' = + 12), to' = |(rii — "2)- Since the rotation with respect to J~ a conserves the quantum number j as 
noted in eq. Q, we also find that 

k + s = 2j ; I + s = 2f . (16) 

Taking into account that to = —j, — j + 1, • • • , j and m' = — j', — j 1 + 1, • • • , j' for a given j, j', respectively, it follows 
that eq. lfl5|l can be reexpressed as 



P« (*) = E E E c ^ • 4,1*' • e-*C B '+^-«- B i'+»'.i'-»')t • {j^-jle-^/^rMU'irn'l^^lj'-J-f) 

j,j'=0 ro=— j m'=—j> 

(17) 

where fe = 0, 1, • • • , 2 j and Z = 0, 1, • • • , 2j', respectively. Then, in defining k = fc — j and 1 = Z — j we relabel p$ 

in terms of /3^j , where k = —j, + ,j and 1 = — j', —j' + 1, ■ • • , j' , respectively. This clearly shows that two 
coupled harmonic oscillators can naturally be treated as an angular momentum oscillator. We now make use of the 
coupled-boson algebra to evaluate the reduced density matrix p^(t) below for two physically interesting, transparent, 
and mathematically manageable initial states, |^(0)), to illustrate the entanglement measure M^t)) 1 namely, the 
case of a simple product state, \Ni) | AT 2 ) (case I), and the case of an eigenstate of the exact Hamiltonian in eq. (J3j) , 
\Ni,N 2 ) H of eq. ifTTj) (case II). In terms of physical interest, the simple product state in case I is an initial condition 
which inherently assumes no entanglement at time zero, whereas case II, using the exact eigenstate of the Hamiltonian, 
implicitly includes entanglement. 

Let us begin with the case I: \ip(0)} = \Nx) I-ZV2). In the coupled-boson representation, this state becomes |J;M) 
where J = ±(N 1 +N 2 ) and M = \{N l -N 2 ). We now have c j>m = (J; M\ j; m)* H -8 3 j and c*, m , = (J; M| j'; m!) H -b jl j 
in eq. 1|17|1 which then reduces to the diagonal form, 

J 

/£?(*) = <5ki E ( J 5 M I J ' m )H (J;M\ J\m') H x <J;k| J;m) H (J; 1| J;m% e^ m ' m '^\ (18) 

m,m' =— J 

where 5\a results from eq. lfl6ll with j = j' = J. In eq. (fl8]) . all the matrix elements can be written in terms of 
(J;/i2| J;fJ,i) H = ^1/^2(7) with (fix — m,m!) and (/i2 = M, k, 1), where UjiP^fr) is formally explicitly expressed as^ 



UP M / (JT^j! ( sin i) Ml ^(cosj) ^ (^L\ J -^ /fcos-y+iy+^fcos-y-iy-^) 

- ( 2J+ , 2)v / (J + /il)!(J _ /il)! Ucos 7 J {(C0S7+I) (cos 7 1) ) 



(19) 



here we see that ^2^1(7) — ujti, U ( — 7) > an d {^1/^2(7)}* = ^ii*a(7)< Particularly for the case that /i2 = J, this 
reduces further to 



which will be used later. Equivalently, U^P^ (7) can be expressed in terms of the tabulated Jacobi polynomial^ 
V^'^\x) in the form 



rj(J) 



'HM = (-1)*+* y gggpg ( sin 2)— x (oob2) w+w • {^-^^'(oobt)} • (21) 
Reexpressing pj^' in eq. lfT8|) in terms of C^u.i/h 3 (7)j we obtain 

pLV(*) = Ad E {^iW}' • ^m(7) • ^(7) • {^1(7)}* • e^"'^ = <5 kl • /«(<) • {/«(*)}* = £ k . • |/^(t)| a , 



m,m' — — J 

where 



(22) 

= £ ^/(7)-^(7)- e - imflt (23) 



m— — J 



5 



with the property that f^{0) = 5m1c [here, we used the relation that {f^2f(7)}* = f^2f(7)]j f$p{t) can be 
explicitly evaluated for any set of the harmonic oscillator initial condition numbers, {Ni, N 2 }. We refer the reader to 
the Appendix [A] for the explicit expression of eq. (|23j) obtained from the substitution of eq. (|2l)l . which will be useful 
below. From (|22| we obtain the measure of entanglement given by 

J 

M im) = 1- E l/k } (*)l 4 > °- ( 24 ) 

k=-J 

We note generally that in eq. |23|) is a time-dependent oscillatory function with a modal distribution coefficient 

^m2f M ' ^mk(7) governing the oscillatory strength. This clearly influences the measure of entanglement as noted in 
eq. (|24j) and as seen in Fig. Q]for specific values of J. 

Let us now apply eqs. f22|) and f2~4|) for specific values of J. For J = h, the total occupation number A^i + A^ = 1, 
there are two possible simple product states | \\ |) , |^; — 5) describing this initial condition. We can easily evaluate 
p (1) (t): for |-0i(O)) = ||; |), from eqs. ([20]) and J23]) we obtain /^(t) = cos ^ + 1 (COS7) sin ^, and then the 2 x 2 
diagonal matrix from eq. l[22|) . 

pH _ t = I (sin 7 ) 2 {1 - cos (to)} ; = 1 - ^> . (25) 

2 ' 2 2 ' 2 2-2 

Similarly, for 1^2 (0)) = |§; — |), we acquire the diagonal matrix p^\t) whose elements are exchanged from eq. 

(|25)l . Due to the fact that ^[/^^(i)] = P^ 1 ^*)] < 1, the reduced density matrices, Pi (t) and p 2 ^(i), represent 
mixed states, respectively, with both directly reflecting the appearance of entanglement between the two coupled 
oscillators; this is noted in Fig. Q] with the properties of M^u)) an d J^i\^ 2 (t)) as given in eq. 112411 . This appearance 
of entanglement in the time evolution from the disentangled initial product state for each case is clearly attributed to 
the interaction Hamiltonian in (!]), i.e., Hk (a\ a 2 + Si a 2 ). Likewise, for J — 1, |, 5 (or iVi + = 2, 3, 10), we also 
obtain the entanglement in the time evolution resulting from the interaction Hamiltonian. Due to their complex forms 
of eqs. ([22]) and lf24j) [see also (|Al[) and (|A2|) ]. respectively, we simply plot the exact numerical results of -M\^(t)) f° r 
N = Ni + N 2 with |-0(O)) = |iV> jO) = | J; J) (see Fig. [J). We clearly see that M\^ t )) increases with N (or J). 

Two points deserve comment here. First, we see that for a given total occupation number A~ = 2J [constant of 
motion; cf. Q], M-\Mt)\ is periodic (see Fig. Q]), as is the reduced density matrix p^'(t) due to the clear oscillator 

dependence on time as seen in in ((23]) ; on the other hand, this is not the case for a single oscillator coupled with 
a thermal reservoir modeled by a sea of infinite oscillators^ 2 -, for then, dissipative decay of its occupation number into 
the the reservoir is effective. Second, in calculating p^ (t) = Tri p(t) and comparing results to the corresponding 

expressions in (fl4|) - (jTHJ) , one easily arrives at the fact that pff = fP^ _, for any J, which confirms the relationship 

valid for any system state \ip{t)) that V[p^} = V[p^], thus indicating the equality of mixture. 

Next, we consider case II, in which the initial state is given by an eigenstate of the exact Hamiltonian in eq. ([3]), 
\Ni,N 2 ) H = \ J;M) H with J = i(ATi + N 2 ) and M = \{N X - N 2 ). From eq. (fT7j), we easily obtain the reduced 
density matrix and the measure of entanglement, 



= ftd • (.7; M e W» J; = Su ■ |/^T , (26) 
where 

/f = ^(7). (27) 

Here, is time-independent, whereas ffp(t) from eq. ([23I is not as obtained in the previous case. Also, we can 
easily show that p^ = /3_ k _ v Then, as in eq. I|24p . we have 



2 



f(2)|2 



M m)) = 1 - E l^k I 4 - (28) 



k=-J 



Particularly for M = J, namely N 2 = 0, from eq. ([201 with UjQ(j) = Ufcj(—j), we easily find that 



/f = ^(7) = VmtP^T (cos*)* (s^) 2J - k , (29) 
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where k — J + k. Therefore, from eq. (|26|) we arrive at the fact that the occupation number k in oscillator 1 is 
described by the binomial distribution B(k; 2J,Pi) with a trial probability p% — (cos^) 2 , which remains unchanged 
in the time evolution. Along the same line, for oscillator 2, we also have 

/f - <(7) = (sini)^ (co si ) 2J - fc (30) 

and, from eq. |26|) . the binomial distribution B(k; 2J,p 2 ) with p 2 = (sin ^) 2 for the occupation number k in oscillator 
2. These reduced density matrices, and j^ 2 \ represent mixed states for 7^0, respectively, due to the fact that 
^[Pi^] = ~P\fi2 ] — 1) thus indicating the entanglement between the two linear oscillators, < M.\^it)) obtained from 
eq. (|28f . Furthermore, as in the previous case, M\^[ t )) increases with J (see Fig. [21); it is noted in the figure that 

■M-Wt)) — ► 1 as J — > oo. 



3. INTERACTION BETWEEN A LINEAR AND A NON-LINEAR OSCILLATOR 
3.1. Quantum-classical behaviors in the Heisenberg picture 

The second coupled system under investigation consists of a linear oscillator and a (non-linear) angular momentum 
oscillator. This system approximates the interaction between a field mode and an atomic TV-level system under ide- 
alized conditions resulting from neglecting dissipation and the coupling between other surrounding atomic systems. 
Following the seminal work of Ref. i, which considered the analysis quantum- mechanically and also classically to 
compare one case with the other systematically, it is convenient to employ, for the linear oscillator, dimensionless 
coordinate x and momentum p, where x = x and p = P , satisfying [x,p] p = i, and, for the angular mo- 

mentum oscillator, dimensionless angular momentum variables j x , j y , j z with [j r , j s ] = ie rs tjt, where the brackets 
[, ] p represent the commutators for the quantum-mechanical description, and i times the Poisson brackets for the 
classical description, respectively. 

The Hamiltonians of the two individual oscillators are then expressed as 

Hi = ifiwi (x 2 + pj) , H 2 = hu 2 (j,) 2 , (31) 
respectively, and the interaction Hamiltonian is given in the rotating wave approximation by 

H12 = Uk (a x (j + ) 2 + a\ (j_) 2 ) , (32) 

leading to the total Hamiltonian 

H = Hq + H\2 ; Ho = Hi + H 2 . (33) 

Here, we have 

ai = 75 ( x i +*Pi) . a \ = 72 ( Xl ~ * Pl ^ ' = 7f {(j^)2 + *(jy) 2 } . (j-)2 = 73 {(ix) 2 -i{jy) 2 } ( 34 ) 

(classically, Hermitian conjugation corresponds to complex conjugation). These non-Hermitian variables and (j z ) 2 
clearly obey 

= 1, [(j+) 2 , (j-) 2 ] P = (J*)a - [(J+) 2 . (j») a ] p = -(j+) 2 . [(j-) 2 . (j*) 2 ] P = CJ-)a ( 35 ) 

with all other brackets vanishing. From now on, let us restrict ourselves to the resonant case, wi — uj 2 = u. In order 
to make later calculations simpler, we then introduce the reduced variables^ specified, with the aid of the identity 
((5]) [for the classical description, the commutators therein have to be obviously replaced by i times Poisson brackets] 
and the rules in <f35|l . by 

Ai = e~ iH ° t l h ax e lH " t / h = ax e luJt , A\ = er lH ° l l h a\ e tH ° l / h = a\ e~ luJt (36) 

and, similarly, 



ax, a[ 



(J+h = (j+)a e"^* , (J-) 2 = 



(j_) a , {J z ) 2 = (j,) a , 



(37) 



7 



respectively. These reduced variables are easily seen to satisfy the same bracket relations as the corresponding 
unreduced variables, namely, 



A U A\ 



i, [(j+) 2 , 



(J,) a , [(J+) 2 , (J*) 2 ] p 



(J + ) 2 , [(J_) a , (J z ) 2 ] p 



(•/-) a 



The interaction Hamiltonian in (|32|) is now rewritten as 

tf 12 = fi/S {Ai (J+) 2 + 4 (J_) 2 } , 

and the equations of motion in terms of the reduced variables will be given by 

ihO(t) = [0(t), H 12 ] p , 

where O stands for any of the above reduced variables. 
From eqs. p8 |l . pQ" |> it follows that 



(38) 



(39) 



(40) 



Ax = —in (J-) 2 ,4 — in (J+) 2 j 
(J+) 2 = -i/t4 (J z ) 2 , (J_) a =«4i (J z ) 2 , (J 2 ) 2 = -in (J+) 2 - 4 (J_) 2 } 



(41) 



and then 



(J z ) 2 = -n 2 {(J_) 2 (J+) 2 + (J+) 2 (J_) 2 + (Ax4 + A\A,) (J 2 ) 2 } = -2 k 2 {(J_) 2 (J+) 2 + Ax4 (J,) a }, (42) 

which, as shown previously,— hold both quantum-mechanically in the Heisenberg picture and classically on the basis 
of the Poisson bracket in phase space. 

From the equation of motion in eq. f40|) it turns out that the sum of the (dimensionless) energy of both oscillators, 
E = ni + (Jz) 2 w ith n i = 4 ^i (without the zero-point energy of the linear oscillator in the quantum-mechanical 
description), and the interaction Hamiltonian H12 are constants of motion, respectively, which are determined by 
a given initial state. As well, the square of a given total angular momentum of the angular momentum oscillator, 
J 2 = (^01)2 + (Jy) 2 + (^2)2 = (^-)2 + (J-) 2 (J+) 2 + (Jz) 2 > is also a constant of motion. Then, substituting the 

expressions 



A X A\ = E - (J,) a + [A u 4] , 2 (J_) 2 (J+) 2 = J 2 (J z )' - [(J+) 2 , (J_) 3 ] 
(note that there are no subscripts p of the brackets) into eq. {42]), we obtain 

(J z ) 2 = k 2 {3 (J z ) 2 2 - 2 (E + [A,, 4]) (J,) 2 + [(J+) 2 , (J_) 2 ] - J 2 } . 



Let us introduce the notation 



Ai,A\ 



A! , [(J+) 2 , (J_) 2 ] - A 2 (J z ) 2 



(43) 



(44) 



(45) 



where A* = 1,0 with k = 1,2 correspond to the quantum-mechanical and the classical description of each oscillator, 
respectively. Eq. f44|) is now reduced to 



(J z ) 2 = k 2 {3 (J z f 2 - (2E + 2 Ax - A 2 ) {J z ) 2 - j (j + A 2 )} , 



(46) 



where the constant of motion J 2 = j (j + A 2 ) with j being classically (A 2 = 0) the total (continuous) angular 
momentum and quantum-mechanically (A 2 = 1) the corresponding quantum number. Here, {J z ) 2 is, clearly, an 
operator for A 2 = 1, with the property ((J z ) 2 ) ^ {(J z ) 2 ) 2 in general. Eq. (J46j) is a non-linear differential equation 
for {J z ) 2 in both the classical and quantum-mechanical description. Since [(J z ) 2 ,E] p = [(J z ) 2 ,E] p = (note that 
E = n\ + (J z ) 2 ), it follows that in the quantum-mechanical consideration of eq. f46|) the constant of motion E can be 
treated as a c-number. 

We now consider eq. (f46]) with two interesting initial states. First, for the ground state with n\ = and (J z ) 2 = —j, 
then E = E 9 = — j; it immediately follows that eq. l(46|) . at t — 0, becomes 



2 k 2 (Ai - A 2 ) j 



(47) 
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Quantum-mechanically, this state corresponds obviously to \tp g ) = |0) 1 \j, — j) 2 . From the fact that each oscillator 
with the minimum energy cannot give up further energy, we necessarily have (J z ) 2 — here, which indicates that 
only Ai = A2 = 0, 1 simultaneously are physically allowed. For the quantum description, the zero-point fluctuation 
Ai = 1 would yield a positive force on the angular momentum oscillator, and the fluctuation A2 = 1 of the angular 
momentum oscillator would lead to energy transfer into the linear oscillator; however, the two formal processes are 
always cancelled such that they are not part of the real (or measurable) physical processes of energy transfer.— 
Therefore, the "semiclassical" option with Ai = and A2 = 1, and vice versa, is not physically admissible for the 
ground state \ip g ), and is thus inappropriate for the system analysis. 

Another initial state under investigation is the state E = j with n\ — and (J z ) 2 = i> quantum-mechanically 
I Vj ) = |0)i I J , J? ) 2 ■ Then, eq. |46|) . at t = 0, reduces to 



t=0 



(48) 



independent of A2! This demonstrates that only for the quantum description of the linear oscillator (Ai = 1) is 
spontaneous emission available, which arises from the zero-point fluctuation. Interestingly, it points out that for 
Ai = 0, there exists an unstable equilibrium leading to the absence of spontaneous emission; the unstable equilibrium in 
the classical description (Ai = A 2 = 0) and the presence of spontaneous emission for Ai = 1 were thoroughly discussed 
in Ref. d. Here, we connect the spontaneous emission process with the notion of entanglement. In considering the time 
evolution of the wave-function for the spontaneous emission in the Schrodinger picture, |?/>(t)) = e~ lHt l h \j,j) 2 
with the Hamiltonian H in eq. (|33)l . we now necessarily have the entangled state for t — > + , i.e., early time beyond 
the initial time, which is given by 



\m) = |0>i \j, i) a - it {« (j + |) 10)! \j,j) 2 - K VJ|1>! \j,3 - 1> 2 } + 0(t 2 ) . 



(49) 



This entanglement is explicitly shown in a compact form for the two-level case (j = ^) later in eq. (|58|) and in the 
more complex form for the three-level case (j = 1) in (|7Tj) . Therefore, for the quantum description Ai = A2 = 1, 
the entanglement is always temporally present in the spontaneous emission process. Here, interestingly enough, we 
still have the spontaneous emission for large enough j, which immediately leads to the fact that the entanglement 
persists even as the coupled angular momentum oscillator is taken to the limit of a large number of levels, a limit 
which would go over to the classical limit for an uncoupled angular momentum oscillator. In essence, the state \ipj) 
of the uppermost excited state of a j-level angular momentum oscillator is always coupled to the vacuum state of the 
linear oscillator. 

From eqs. 1|39|1 . (J4TJ) we also obtain 



(J z )l + K 2 = 2n 2 {a x A\ (J+) 2 (J_) a + A\ A, (J_) a (J+) 2 } , 
where K = Hi 2 /h, and subsequently, with the aid of eqs. I|43p . (|45| . the expression 

(J z )l = 2k 2 [{(J z ) 2 - (E + %)} {{J z )l - + A 2 ) - |AxA 2 )}] + (E + 4^) Ai A2 k 



K 2 



(50) 



(51) 



In Ref. [1, it was indicated that in the classical description (Ai = A 2 = 0) eq. (HJ) has the form of an equation for 
the vertical position J z of a classical spherical pendulum, where the constant of motion K represents an angular 
momentum about the vertical axis through the center. It turns out that for m = and ( J z ) 2 — j we get K = from 
(j39"l) with Ax = y/n{~e iei , A\ = ^fn{e- l(>1 and J± = ^(J 2 - (J z ) 2 )/2 e ±lf>2 (note that in the quantum description 



the constant of motion K 2 7^ 0) and then (J z ) 2 = from f5lj) . With (J z ) 2 = in eq. J48|, this reveals that 
the initial state (ni — 0, ( J z ) 2 = j) corresponds to an unstable equilibrium of the pendulum. 

Here, it is also worthwhile pointing out here that the semiclassical treatment of Jaynes-CummingsAS reproduces this 
spontaneous emission for j = 2 only, with the same decay rate as that given by the quantum-mechanical description; 
this is accomplished by adding a phenomenological damping term to the classical equation of motion of the linear 
oscillator which would result from the interaction with the instantaneous expectation value of the dipole moment of 
the angular-momentum oscillator. Therefore, the Jaynes-Cummings model just provides an "artificial" picture of the 
spontaneous emission without describing the actual physical processes involved therein (e.g., entanglement), and also 
without offering a direct way for its extension to more than the two- level case (j > 1). 
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3.2. Appearance of the entanglement in the Schrodinger picture 

In order to study the appearance of the entanglement between the linear and the angular momentum oscillators in 
the quantum description, we discuss the state evolution of the total system in the Schrodinger picture given by 

H = hw{n l + |) + huj J z + Hk(A 1 J + + A\ J_) (52) 

from eqs. f31j) - ([331 with the reduced variables in eqs. (|36| - l(37|) for the resonant case, u>\ = u>2 = w. We will 
explicitly consider below the measure of entanglement for the cases of j = |, 1, | of the angular momentum oscillator. 
To this end, as will be seen, the (2j + 1) x (2j + 1) irreducible matrix representation for the total system is employed, 
which is also systematically applicable for j > 2. Here, we assume that the given total energy (a constant of motion) 
is given by E = n% + J z = ni + j > j. 



3.2.1. Casej = \ 

We first consider the case of j = 4, namely, a two-level angular momentum oscillator interacting with a linear 
oscillator. For the total energy E = n\ + j = n\ + \, the basis {\m + 1, —5), \ni, \)} allows us to rewrite eq. (|52|) in 
explicit operator form 1 ^ as 

h = {m + \)hu { |m + 1,-5) (ni + 1, -|| + \m, \) (ni, || } + 

+ 1 { |ni + l,-|>(m,i| + |ni,|>(ni + l,-||} . (53) 
This can be written in a 2 x 2 matrix form as 

H = (m+ 1)W1 2 + -j=hKy/m + la x , (54) 

where a x = 2J X denotes the Pauli matrix. Using the identity, 

e la ^ a x e- <a *" = (cos2a)o- 2; + (sin2a)o- z (55) 
with a — j , we arrive at the diagonalized Hamiltonian for eq. |54|) as 

H d = (ni + l)Hu>l 2 + ^y / n 1 + l(sin2a)a z . (56) 
Accordingly, we acquire the energy eigenstates as 
|1 2 ) = e -<f*. | ni + i,_i) = _^(|„ 1 + i s _i) + | ni> i)) , |2 2 ) = e -<*»» \n x ,\) = ^(-K + l,-±) + |m, §» (57) 

with the corresponding eigenvalues, E\ = (ni + 1) So; — ^|-\/ni + 1 and £2 = [m + 1) ftw + ^= v / n~H _ X respectively. 

Then, just as in discussion related to eq. (fl3| for the two coupled linear oscillators, the state evolution for the 
initial state | -0(0)} = Ini, |) is given by 

m)) = 71 11 2 > e-* Blt - ^5 |2 2 ) e-t*" . (58) 

This state explicitly includes the possibility of spontaneous emission (rii = 0). From eq. (|58| . we easily acquire the 
reduced density matrix of the angular momentum oscillator as 

p^{t) = Tn m))W)\ = sin2 (^V^7TT) l-i) (-I) + cos 2 (^V^TTl) ||> (I) , (59) 

which represents a mixed state. Therefore, it immediately follows that 

V[p^] = sin 4 (^V^TTT) + cos 4 (^V^TTT) , (60) 

with M.\Mt)) = 1 — V[p^] which oscillates between and \ in the time evolution, thus indicating the appearance of 
entanglement between the linear and the angular momentum oscillators (see Fig. [3]). Furthermore, from eq. lf59"l) we 
find, after a minor calculation, that 

{Jz) = Tr 2 { P W(t) ■ J z } = \ , cos(^^7TT) , (J 2 Z ) = \ , (J 3 Z ) = I cos(^|^7TT) • (61) 

Therefore, from eqs. l(46|) . l(6T|) . it follows that 

(J z ) = -2( ni + l) K 2 (J z ), (62) 
which obviously yields the consistent result with that in eq. (|48|) for the case that n± = and J z = j = |. 
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3.2.2. Case j = 1 



Next, we consider the case of j = 1. For the total energy E = n\ + 1, from eq. Il52[l with the basis {|ni + 2, —1), \n\ - 
1, 0), \n\, 1)}, the Hamiltonian takes up the 3x3 irreducible matrix representation 



H 



/ (ni + §) hw V«i + 2 hK \ 

\Jrt\ + 2 Hk (ni + |) hw yjn\ + 1 Hk 
\ ■ s /n 1 + THk (ni + |) fiw / 



(63) 



Based on the generating operators of the group 5?/ (n) ,— where n = 2j + 1 (see Appendix [B]) , the Hamiltonian in eq. 
(j63"l) can be reexpressed as 



H = hui (ni + 



Is 



H12 , H\2 = thKy/m + 2ui2 + hn \Jri\ + 1^23 



(64) 



where U12 and U23 are generating operators of SU(n), and their matrix forms are explicitly given in eq. l|B3j) , 

Here, it is interesting to note that for a large «i > 1 such that ^/m + 1 « ^/ n i + 2, the Hamiltonian in eq. ([62 
can be transformed approximately to 



H 



(ni + |) fiwl 3 + ^2 (ni + 2)hnJ x . 



(65) 



which formally corresponds to the exact Hamiltonian for j = \ in eq. (|54| with j x = (T 2; /2. Thus, for a large m, we 
have an approximate diagonalized representation Hd for j = 1, based on the previous result. 

For the case of general n±, proceeding from eq. (|63|) . we easily obtain the energy eigenstates with respect to the 
basis {|m + 2,-1), |m + 1,0), |ni, 1)} as 



/ v«i + 2 \ 



113) 



y/2 (2m + 3) 



V2ni + 3 



-1 



V2"i + 3 



/ a/th + 1 




, |3a) = 



\ -V^T+2 



^/2(2m +3) 



/ V^T2\ 

-V2m + 3 



and the corresponding energy eigenvalues, 

Ei = (m + |) huo + Kh^2n x +3 , £ 2 = (ni + |) Kw , £3 = (m + I) hw - Khy/2nTT3. 
respectively. Therefore, the diagonalized Hamiltonian is given by 

H d = (ni + |) 1 3 + ft K V2m+3(|l 3 )(l 3 | - |3 3 )(3 3 |) . 



(66) 



(67) 



(68) 



In order to diagonalize the Hamiltonian H in eq. lf52"|) for higher manifolds of j, it is instructive to systematically 
analyze the above diagonalization of H and confirm for j = |, 1, In seeking the rotational operation that diagonalizes 
H in eq. lj64"|) to achieve Hd, it is necessary to consider the unitary, Euler-like transformation 



_ g2«3 W 13 g«C«2 «23 g«ai Di 2 jj g — IQl t>12 i 



(69) 



comprising three sequential, non-commuting rotations with respect to (£12, #23, #13) defined in (jB4|) . with Euler angles 
(ai, «2) ai 3 ). Then, as previously accomplished for the case of j = 5 in eqs. (|54|) - (J56j) , here we similarly choose 
(ai 



-f;cosa 2 = + l/v2ni + 3, sina 2 = — v^i + 2/y/2n\ + 3 ; a 3 = — |) to obtain iJ^ in eq. (|68|) . 
This clearly shows that the diagonalization of the Hamiltonian in eq. lj64"|) cannot be described by a simple rotation 
such as i? (j=1 " ) (a, 0, 7) = e~ laJz e~ l ^ Jy e~ llJz of an angular momentum oscillator^ whereas for j = h, the rotation 
i?W = s)(0, |, 0) was used in eq. (|55|) to diagonalize 77 in eq. (J54J) . However, for the case of j = 1 and large ni > 1 
given in eq. lj65|) . the diagonalization can be approximately accomplished by i?W (0, f,0). Finally, it can be shown 
by induction, for general j, that the diagonalization of the Hamiltonian in eq. lf52|) is generally characterized by 
\(v? — n) rotation angles {a\, a<i, ■•■ , ai( n 2_ n )} with respect to the set {vjk} as given in eq. (|B2[) . where n = 2j + 1. 

Let us now consider the state evolution with the initial condition |^(0)) = \rii, 1); just as in eq. ([13|. from eqs. 
(I66l). (1671) we find that 



W)) 



2(2ni+3) 



|1 3 ) e-n 



Bit 



ni+2 
2ni+3 



|2 3 ) e" 



-E 2 t 



ni+1 
2(2ni+3) 



|3 3 ) 



E 3 t 



(70) 



11 



From this state, we obtain, after some calculations, the reduced density matrix of the angular momentum oscillator, 



P (2 \t) = TnlVWXVWI 



(2ni+3) 2 



4(rn + l)(m + 2) sin 4 (^±±2 K ^ | — 1><— 1] -f- 



(ni + l)(2m + 3) sin 2 (y/2m + 3 nt) |0> <0| + |2(m + 1) cos 2 ( V2 " 2 1+3 Kt\ + l} 2 |1)(1| 



(71) 



This reduced density matrix represents a mixed state with P[p( 2 '] < 1 and M.\Mt)) > 0, directly reflecting the 



entanglement of the total state \ip(t)) [here, eq. ([7Tj) for j = 1 is comparable to eq. (|59|) for j 



Furthermore, by 



noting that |^>(i+T)) = e 



\tp(t)) = \ip(t)) for a given t with the diagonalized form Hd in eq. (f68|) . we can easily 



show that the state evolution in eq. l(70|) displays periodicity with period T 



2tt 



as noted in Fig. [3l This means 



V2ni+3 k 

that the time evolution operator U(t) = e - lHt / h ^ expressed here with respect to the basis {I3, \ p \p = 1, 2, • • • , 3 2 — 1} 
given in Appendix [Bj as 



u(t) = i 3 + u p& v 
P =i 



(72) 



where U p (t) = Tr {U(t) A p } and U p (0) = for all p, has the periodic property that U p (t) = U p (t + T). 



3.2.3. Case j = § 

Let us consider the case of j = |, where the Hamiltonian of eq. (|52| for the total energy E = n\ + § is given by 

i? = hcu {n x + 2) i 4 + £12 ; H 12 = h K (^J^^u 12 + ^2( ni + 2)u 23 + ^Mnitli^ . (73 ) 

The diagonalized form of eq. ([731 can be obtained, after tedious calculations similar to those described in the previous 

H d = huj (m + 2) i 4 + f { VeTTe; ( |1 4 ) (1 4 | - |2 4 ) <2 4 | ) + a/^T - ^ ( |3 4 ) (3 4 | - |4 4 ) (4 4 | ) } , (74) 

where Si = 20 + 10ni; 9 2 = 2^73 + 64ni + 16n 2 , and the |p 4 )'s with p = 1, 2, 3, 4 denote energy eigenvectors. 

Due to the complex form of the reduced density matrix of the angular momentum oscillator, p^ 2 \t) — 
Tti\^(t)){ip{t)\, and the measure of entanglement, M\^ t )) = 1 - V[p (2) {t)}, obtained from \i/j(t)) = e - %Ht / h |^(0)) 
with IV'(O)) — \ n i: |)j we simply plot the exact numerical results of -M-UjiCt)) in Figs.[H[5j 

Based on the numerical analyses^ and the statistical approximation in terms of classical random variables^, it has 
been shown that (J z (t)} with |^>(0)) = |0, j) exhibits aperiodic behavior in its time evolution for all j > |, even for j 
large enough (i.e., in the classical limit for an uncoupled angular momentum oscillator), while the classical counterpart 
displays periodic motion on the spherical pendulum [cf. eq. ([5Tj) with Ai = X 2 = and K = 0; note that for this initial 
state, the angular momentum K about the vertical axis through the center vanishes]. This shows that the aperiodicity 
in the time evolution would be of non-classical origin. We indicate this aperiodicity for any n = 2 j + 1 > 4 here by 
noting the impossibility of having the periodic property U(t) = e ~ lHt / h — t n for arbitrary t; from eq. ([74[l . we find, 
after some calculations, that 

U(t) = e~i" dt (75) 
= e -«rf(m+2) ( e -i§ Vd+^t |i 4 )(i 4 | + e lJ i VBl+ ° 2t/2 \2 4 ){2 4 \ + e-^ V0l - &2t \3 4 )(3 4 \ + e'* V«i-»»t |4 4 )(4 4 | ) . 

Since here the phase factors, n \fW\ + 9 2 /2 and K\/Wi — 9 2 /2, are incommensurable with respect to each other, we 
immediately note that not all Uk(t) in the expression 

4 2 -l 

U(t) = U + J2 U k(t) A fc (76) 
fc=i 

can simultaneously vanish periodically, which leads to the aperiodicity in the state evolution \i/)(t)) = U{t) |"0(O)) (see 
Figs. [H [5]). Along the same line, it can be shown that we have the diagonalized Hamiltonian for any j in form 

H d = hoj (ni + i+j) l n + Kh{/3i (|1„)(1„| - |2„)(2„|) + ••• + 

(3 m ( I (2m - 1)J ((2m - 1)„| - |(2m)„) ((2m) n | )+•••+ /3 n/2 ( |(n - 1)„) ((n - 1)J - |n n ) (n n | ) } ,(77) 
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where n — 2j + 1, and (3 m — (3 m {n\) with m = 1, 2, • • • , ^; |p„)'s denote energy eigenvectors. Here, each (3 mi is, in 
general, incommensurable with any other (3 m2 , where mi ^ m 2 , and for the case that n is odd, one of /3 m 's always 
vanishes. Therefore, we have 

U(t) = e -**<* = e- Mt( " 1+ ^ +j) (e-' iK/?lt |l„)(l„| + e lK ^ t \2 n )(2 n \ + ••• + e~ iK ^ * |(2m - l) n )((2m - l) n | + 

e «A»*|(2m) n )((2m) n | + • ■ ■ + e"™ <W |(n - l)„)((n - 1)J + e™"-/"* |n„)(n„| ) , (78) 

thus showing that U(t) ^ 1„ for any t unless all f3 m are commensurable with each other. From this, we easily find 
that the aperiodic behavior in the state evolution, \ip(t)) = U(t) |-0(O)) , will survive and even increase with n > 4, i.e., 
j > | [note that for j = 1, we obviously have (3\ only, which leads to the periodicity in the state evolution: U{t) = I3 
for the case that k(3\ t = ±2tt, ±47r, ±677, • • •], 

4. CONCLUSIONS 

In summary, we have investigated the fundamental dynamics of two interacting oscillators: in one scenario, two 
linear oscillators, and in the other scenario, a linear and a non-linear oscillator have been considered. For the first 
scenario, based on the coupled boson representation we compactly described the equation of motion in the Heisenberg 
picture and also systematically studied the quantum entanglement in the time evolution of the total wave-function 
\ip(t)} for various initial states in the Schrodinger picture. From this, we found that the appearance of entanglement 
in the time evolution from the disentangled initial simple product state is attributed to the interaction between the 
two individual oscillators. Also, the measure of entanglement increases with the total occupation number. 

For the second scenario, quantum versus classical behaviors have been studied, based on the Heisenberg equation 
developed in terms of relevant reduced kinematics operator variables and parameterized commutator relations. By 
setting the corresponding commutator relations to one or zero, respectively, the Heisenberg equations are shown to 
describe the full quantum or classical motion of the interaction system, thus allowing us to discern the differences 
between the fully quantum and fully classical dynamical picture. In addition, for this second scenario, in the fully 
quantum- mechanical description, we considered special examples of j = ^,1,§ for the coupled angular momentum 
state, demonstrating the explicit appearances of entanglement. The entanglement increases with j and so persists even 
as j — + 00, a limit which would go over to the classical picture for an uncoupled angular momentum oscillator. This 
entanglement occurs because the uppermost excited state of the j-level angular momentum oscillator is always coupled 
to the vacuum state of the linear oscillator, a purely quantum coupling which manifestly gives rise to spontaneous 
emission. We have also shown that the dynamics of this scenario can be systematically described by the irreducible 
matrix representation based on the generating operators of the group SU(n), while that of the first scenario can 
be given, based on the coupled boson representation, simply by the rotations of an angular momentum oscillator. 
For the coupled linear-angular momentum oscillator, this system was shown to display periodicity in the measure of 
entanglement for 3 — \ and j = 1, whereas for j = | and beyond, the measure of entanglement was shown to be 
aperiodic; this aperiodicity is apparent from the form of the diagonalized multi-level Hamiltonian and the resulting 
structure of the time evolution operator. 
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APPENDIX A: MATHEMATICAL SUPPLEMENTS TO EQ. ([23]) 
From eqs. (HJ) and {23]) with the relation that ^^(7) = ^km( - 7)> we find that 

/^w = (-i) M v /I M{ii r ( sin i) A " k ( c -i) M+k >< 

J2 (-l) m (cos|) 2m • {^M m,M+m) (cos7)} • {^ m m k '" i+k) (cos 7 )} • e-"" flt . (Al) 
m=— J 
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Particularly for the case that 7 = §, namely u>i — u>2, eq. (|A1[) reduces to 

J J-M J-m 

fl 1] (t) - ir v^SSSf E *~ imQt ■ E E ■ fr) G-M m P ) (V) (/: k J . (A2) 

m——J p—0 q—0 

where (») = 



APPENDIX B: GENERATING OPERATORS OF THE GROUP SU(N) 

The set of the generating operators is given by^ 

A = {\ p \p = 1,2, • • • ,n 2 - 1} = {Ml2,Ul3,M23, ■ • ■ , #12, $13, V23, ,Wl,W 2 ,- ■ ■ ,W n -l} . (Bl) 

where |A| = n 2 — 1, and 



Ukl = Pkl+Plk; i> kl = i(Pkl-Plk); W m = -J m(? ^ +1) (Pll + P22 H \- Pmrn - mP m+1 . m+1 ) . (B2) 

Here, 1 < k < I < n, and \{u k i}\ = \{vm}\ = \(n 2 - n), where {u k i} and {v k i} denote {ui 2 , U13, "23, • • • , «n-i,n} and 
{«12 }, respectively; 1 < m < n — 1, and |{w m }| =n— 1, where {w m } denotes {wi,W2, • ■ ■ ,u>n-i}. 

-Pfei = represents a projection operator for k = I, and a transition operator for k ^ I, For n — 2 (i.e., j = ^), the 

generators (ui2,v±2,wi) exactly correspond to (a x ,a y ,a z )* For n = 3 (i.e., j = 1), three generators {U12, M23, W13}, 
for example, have the matrix representations, 

/ 1 \ / \ / 1 \ 

«ia = Aa + Ai = 1 , u 23 = As + A2 = 1 , u l3 = P 13 + P31 = 0, (B3) 

\oooy \ 1 / \ 1 J 

respectively. Similarly, for 3 generators {1)12,^23,^13}, we have 

/ i \ / 0\ / i\ 

fa = -* , V23 = 01,1)13= 0. (B4) 

\oooy y -i oy \ -i o 0/ 

The set X p \p = 1, 2, • • • , n 2 — 1} can be used as a well-defined basis of the space ofnxn matrices. 
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FIG. 1: The measure of entanglement, y(t) = M\^(t)) given in eq. Il24[| versus time t is depicted for k = 1 and initial states 
1^(0)) = |JV}|0> = | J; J), where N = 1,2,3,10 (i.e., J= |,5) in progressive order from bottom to top; the solid lines are 
obtained for 7 = \ (or Wi = while the dash lines for 7 = \ (or Wi — lu 2 = 2). 
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FIG. 2: The measure of entanglement, y = M.\$(t)), as given in eq. (l28l) , versus J is depicted for /t = 1 and initial states 
1^(0)) ~ \J\ J)h = \N,0)h [cf. eq. I|29p ]: the circles (o) are used to indicate 7 = 5 while the boxes (□) indicate 7 = f . As 
noted in eq. I|27p . y is time-independent. 
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FIG. 3: The measure of entanglement, y(t) = M.\i>(t}) versus time t is depicted for k = 1, Alu = 0, and initial states 
IV'(O)} = |ni, j); the solid lines are obtained for m = 0, while the dash lines for m = 20; j = 1 from bottom to top [cf. eqs. 
lf60]l . iffTjl ]: for j = |, the current yit) takes up the form of the y(t) given in Fig. [TJ for j = 1, the current yit) becomes closer 
in form to the y(t) given in Fig. [l]with the increase of n\. 
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FIG. 4: The measure of entanglement, y(t) = M\^(t)) versus time t is depicted for k = 1, Alu = 0, and \ip(0)) = \ni = 0, j = §); 
we see here the the aperiodic behavior of y(t), and its maximum y m is larger than maxima of the y(t) (solid lines) for j = |, 1 
given in Fig. [3] 
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FIG. 5: The measure of entanglement, y(t) = M\i>(t)) versus time t is depicted for k = 1, Alu — 0, and \ip(0)) = \ni = 20, j = |); 
we also here see the the aperiodic behavior of y(t), and its maximum y m is larger than maxima of the y{t) (dash lines) for 
j = 5, 1 given in Fig. El 



